Creation of a molecular condensate by dynamically melting a Mott-insulator 
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We propose creation of a molecular Bose-Einstein condensate (BEC) by loading an atomic BEC 
into an optical lattice and driving it into a Mott insulator (MI) with exactly two atoms per site. 
Molecules in a MI state are then created under well defined conditions by photoassociation with 
essentially unit efficiency. Finally, the MI is melted and a superfluid state of the molecules is created. 
We study the dynamics of this process and photoassociation of tightly trapped atoms. 

PACS numbers: 03.75.Fi, 42.50.-p, 42.50.Ct 
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The generation of Bose Einstein condensates (BEC) 
of dilute atomic gases has resulted in a remarkable se- 
ries of experiments demonstrating various properties of 
quantum degenerate gases [Q. One of the next major 
goals in this effort is the realization of a molecular BEC. 
A promising route towards a molecular condensate is the 
conversion of an atomic BEC to molecules via photoasso- 
ciation, a process discussed so far for conditions of quasi- 
homogeneous trapping of atomic gases ^ [|, [|, [j| . In this 
Letter we describe a novel path to create condensates of 
composite atomic objects, in particular a molecular BEC, 
based on photoassociation via a Mott insulator state of 
bosonic atoms trapped in an optical lattice ||, Q]. This 
provides an efficient way of generating a molecular BEC, 
avoiding some of the problems encountered in the quasi- 
homogeneous case [0. It also touches upon fundamental 
questions related to the formation of a BEC by "melting" 
of a Mott-insulator (MI) state in a quantum phase tran- 
sition, as opposed to the familiar growth from a thermal 
cloud of atoms ^ . 

Experimental advances in manipulating BECs [[j] , and 
in particular the loading of a BEC into an optical lat- 
tice generated by interfering laser beams have recently 
led to a seminal experiment by I. Bloch and collabora- 
tors [Q. They demonstrated a quantum phase transition 
from a BEC or superfluid (SF) state into a MI by vary- 
ing the lattice laser intensity, as proposed theoretically in 
||. While a SF phase has long range order, the MI phase 
corresponds to the loading of a precise number of atoms 
into each lattice site, i.e. Fock state occupation of lat- 
tice sites. Among the proposed applications of this new 
atomic quantum phase are the study of ultracold con- 
trolled collisions and quantum computing with neutral 
atoms ||. In the present context, the MI phase opens 
the possibility to efficiently create a molecular BEC in 
the following four steps: (i) an atomic BEC is loaded 
into an optical lattice, (ii) the depth Vb of the optical 
lattice is increased to create a MI with exactly two parti- 
cles per lattice site, (hi) a molecular MI state is produced 
by two-color photoassociation of the atoms under tight 
trapping conditions, and (iv) by decreasing the depth of 
the optical lattice the MI state is "melted", and thus a 



molecular BEC is created in a quantum phase transition. 

At the end of step (i) above, we have an ensemble of 
bosonic atoms illuminated by orthogonal, standing wave 
laser fields tuned far from atomic resonance. These laser 
fields generate a potential for atomic motion of the form 
= Sj=i ^Oisin 2 (kxi) with k = 2ir/\ the wave- 
vector of the light and lattice period a = A/2. The 
dynamics of bosonic atoms occupying the lowest Bloch 
band of an optical lattice is well described by the Bose- 
Hubbard model (BHM) |6| which includes the interaction 
U a between particles occupying the same lattice site and 
the tunneling J a of particles from one site to the next. 
The BHM Hamiltonian is given by 

H a = -J a a\aj + eifij + -U a ' s ^hi{n i - 1) (1) 

where is a bosonic destruction operator of a particle 
at site i. The number of particles in site i is given by the 
operator hi = ajdj, and is an energy offset due to an 
external trapping potential. The first sum in Eq. (Q) runs 
over all nearest neighbors denoted by Increasing 
the laser intensity of the trapping laser tends to compress 
atoms near the nodes of the lattice field, and thus leads 
to an increased on-site interaction U a , while the atomic 
tunneling rate J a decreases 0. The BHM predicts a 
quantum phase transition from the SF phase to the MI 
state: according to mean field theory this occurs at the 
critical value Ua ~ 5.8zJ a |l(| with z the number of 
nearest neighbors of each site. This corresponds to a 
(moderate) potential depth of Vb = IOEr where En — 
h 2 k 2 /2m is the recoil energy for atoms with mass m. 

We will first illustrate the dynamics of the BHM Hamil- 
tonian with a time dependent depth Vb(i) of the opti- 
cal lattice controlled by the laser intensity, leading to 
a variation U a (t) and J a (t) in Eq. |l| We assume the 
system initially to be in the SF ground state and calcu- 
late its time evolution for the time dependence shown 
in Fig. |l|a Q. For a model problem of N particles, 
where N is small (sa 10), in a few lattice sites, the time 
dependent Schrodinger equation for the wave function 
<F(i) can be solved exactly. Fig. [I]b plots the eigenval- 
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FIG. 1: a) Time dependence of U a (solid curve) and J a (dot- 
ted curve). The vertical dashed lines separate the SF and 
MI regions expected for adiabatic time evolution of the sys- 
tem in ID. Parameters: we assume 87 Rb and A = 390nm with 
V (t) = V SF + {VMi-VsF)/(l + eyLp{(t 2 -tl)/t 2 s)), t w = 30/ J , 
t s = 40/Jo, Vmi = 20E R , Vsf = 5E R giving J = 0.13E R . 
Jo/z is the hopping matrix element for t — > — oo b) Transi- 
tion from the SF to the MI back to the SF phase for atoms. 
We plot e; (lower two curves are doubly degenerate) against 
t for N — M = 5 in a ID lattice with periodic boundary con- 
ditions, c) Atoms in the SF phase are driven to a MI phase 
(time interval t < 0), converted to a molecular MI phase by a 
Raman pulse (shaded region around t — 0) , and melted to ob- 
tain a molecular BEC (t > 0). We plot e; (lower curve doubly 
degenerate) for the atoms (molecules) before (after) the con- 
version for N/2 = M — 3 in ID, and parameters Jt = J a /2, 
U a = Ub = Uab with time dependence given in Fig. |l|a). d) 
Same as b) but using the Gutzwiller ansatz (the lower curve 
is four fold degenerate). 



ues ei of the one particle density matrix pij = (ajaj) 
(= (i) | aja.y | >]/(£))). As expected, there is one large 

(c) 

eigenvalue for U a < Ua of the order of the number of 
particles. The corresponding wave function is approxi- 
mately given by \ipsv) <x al) N \va.c) for N particles in 
M sites with |vac) the vacuum state. All the other eigen- 
values are small and are associated with the quantum de- 
pletion of the SF state. As U a increases and crosses the 

(c) 

critical point U a all eigenvalues tend towards one, cor- 
responding to a diagonal single particle density operator 
(MI). Upon ramping U a down again the SF is restored. 

An extension of the BHM (Q) in an optical lattice de- 
scribes the situation with atoms and molecules present. 
Denoting the annihilation operator of a molecule by hi, 
we add to the Hamiltonian (|l|) a tunneling Jf, and on-site 
interaction term Ub for molecules, and an atom-molecule 
interaction term of the form H a -b = U a } ) '^2 li b\bia i ai. 
The underlying assumption is that the laser beams gen- 
erate an optical lattice for atoms and molecules with the 
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FIG. 2: Production of ground state molecules using two Ra- 
man transitions. The inset shows the first Raman transition 
(solid lines) and the process where U2 is absorbed before wi 
(dashed lines). For 87 Rb and for a trap frequency of v = 
1MHz, Eq ~ 2 MHz and outer turning point Rtp ~ 500 ao 
(ao=0.0529nm). For the singlet pathway, Eb ~ —32 GHz and 
Rtp ~ 32 ao. 

same structure of nodes and antinodes, although both 
lattices can have different depth. The process of atom- 
molecule conversion by photoassociation is described by 

the Hamiltonian H conv = Cl(t) (b\aiai + h.c.^ /y/2, 

where Q(t) is an effective Rabi frequency which is turned 
on at time t = for a short time interval to convert two 
atoms at one lattice site into one molecule. In practice, 
this process will consist of several Raman steps to go to 
the molecular ground state (see Fig. @). 

We emphasize several distinguishing key features of the 
atom-molecule conversion process in an optical lattice. 
First conversion is most efficient under tight trapping 
conditions (in a regime where tunneling between lattice 
sites is negligible). The high atomic densities associated 
with the strong compression of the atoms opens inelastic 
collision channels that typically quench all lattice sites 
with three or more atoms. Thus we assume that only 
lattice sites with occupation of two (or one) atom sur- 
vive the lattice compression. These loss processes are 
added to our Hubbard dynamics by writing down a mas- 
ter equation for a density matrix that includes these de- 
cay channels. However, since in the MI phase atomic 
number fluctuations are small, there are only a very few 
lattice sites that are actually depleted by this loss pro- 
cess. Second, for (exactly) two atoms trapped at one lat- 
tice site we have a complete microscopic understanding of 
the two-atom dynamics, and the conversion to molecules 
by photoassociation jflj beyond the effective description 
contained in the Hubbard model. 

The two-atom Schrodinger equation |ll| at a given lat- 
tice site separates for harmonic confinement into center- 
of-mass and relative coordinates. The potential curves 
for the relative motion of the two atoms in the trap 
are schematically shown in Fig. 0: for small distances 
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R we have the familiar molecular Born-Oppenhcimer po- 
tentials while the large i?-behavior is dominated by the 
trap confinement. The trapping potential discretizes the 
molecular continuum (scattering) states to form a se- 
ries of harmonic oscillator trap states with frequency 
v = (WqEr) 1 / 2 j% with V the potential depth. The goal 
is to perform a Rabi flop SIT = 7r from the lowest trap 
state of two atoms (i.e. the state associated with the low- 
est atomic Bloch band) to a bound molecular state, with 
O the two-photon Rabi frequency. The discreteness of 
the trap states makes this a bound-bound Raman transi- 
tion. The condition for not exciting any other trap states, 
i.e. to avoid heating, is SI <C v. On the other hand, for 
incoherent processes (as spontaneous decay from the in- 
termediate state) with effective decay rate 7 to be small, 
we must have 7T <C 1. Thus we require a large two- 
photon Rabi frequency and tight trapping, 7 <C SI <C v. 
Note that SI involves a matrix element from the bound 
trap state to a molecular state, which results in a scal- 
ing SI oc ^ 3 / 4 11 , i.e. under conditions of tight trapping 



the two-photon Rabi frequency will be significantly en- 
hanced. We will give specific numbers for these parame- 
ters for the case of 87 Rb below. Thus it is the preparation 
of the two-atom MI phase together with strong confine- 
ment which guarantees the coherent conversion of atoms 
to molecules with essentially unit efficiency. 

Fig. [l]c shows results from the exact integration of 
the Schrodinger equation for three sites with six atoms. 
Starting from the SF phase, the atoms are driven to the 
MI phase and converted into molecules at t — 0. Melting 
of the molecular MI phase then produces a molecular SF. 

To describe the dynamics with a large number of par- 
ticles in 2D and 3D we employ a time dependent mean 
field approximation based on a Gutzwiller ansatz [ [To] . 
For simplicity of writing we consider for the moment 
the case of atoms alone, where we write the wave func- 
tion as the product of superposition states at the vari- 
ous lattice sites, \G(t)) = UZi (£T=o /" ) (*)l"->i)- This 
ansatz is motivated by the success and simplicity of time- 
independent Gutzwiller mean field theory to model the 
ground state and phase diagram of the BHM (l(| . The 
ground state is obtained from the variational principle 
(G\H\G) — fi{G\N\G) — ► min, where fi is a chemical po- 
tential introduced to enforce a given mean particle num- 
ber H . From the time-dependent variational principle, 
{G{t)\ih§- t - H(t)\G(t)) -> min, the following time de- 
pendent equation is readily derived: 
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where = (G\ai\G) = E n /n-iV"/^ is the atomic 
SF density. Eq. (g) is a nonlinear equation for the 
amplitudes fj- l \ which preserves both normalization of 



the wave function and the mean particle number. In 
the SF limit and for a coherent state distribution of 
/W = exp („ |^|2 /2)/V^I, Eq. (g) reduces to a time- 
dependent Gross-Pitaevskii equation for on a lattice. 

By projection to a state with definite parti- 
cle number N, \G N ) = V N \G) / \\V N \G)\\ ~ 

So* di P exp(iA^) ni^i (E5T=0 exp(-irup)f^' |n)i) a more 
consistent description is obtained. The resulting time de- 
pendent Schrodinger equation for the amplitudes /« is 
significantly more complex. It can be shown, however, 
that fi^ of the number projected Gutzwiller wave func- 
tion again obey Eq. (g) , provided the variance of the par- 
ticle numbers at each lattice site satisfies Am 3> 1/ \/~N 
(where An? = (nf) — (rii) 2 with (nf) = J2n n2 \fn\ 2 and 

( n i) — J2n n \fn ^\ 2 )- Note that this excludes the regime 
where we have a precise locking of the particle number, 

i.e., fn = 6 n ,no-> as i R the MI obtained from the non- 
number conserving Gutzwiller for a homogeneous situa- 
tion, when the number of particles is commensurate 
with the lattice sites M and no — N/M. However, a 
precise Fock state is never realized in the time evolu- 
tion we consider since the initial superfluid density is not 
completely destroyed while ramping the optical lattice 
up. Also, an additional trap potential confining the sys- 
tem to a certain region in space ensures the existence 
of a remnant superfluid component. Below we model 
the evolution of an initial SF to an (approximate) MI 
while changing U a and J a by integrating mean field equa- 
tions of the type (g). Fig. |l|d gives the results obtained 
from Gutzwiller theory with initial state given by the 
time-independent Gutzwiller wave function, to be com- 
pared with the exact integration of the time dependent 
Schrodinger equation in ID for a few particles in Fig. gb. 
As expected, mean field theory shows a more pronounced 
phase transition than the few atom ID calculation. 

Using a generalization of the Gutzwiller ansatz to 
include superposition states of atoms and molecules, 

l G W) = 1X=i (Ena,n»=o W|n„,n b )iJ, where n a 
and rib refer to the atomic and molecular occupation, re- 
spectively, we have numerically investigated the creation 
of a molecular BEC in a 2D lattice with a superimposed 
harmonic trapping potential. The results are shown in 
Fig. g. As expected we find that molecules are only cre- 
ated in sites with an atomic occupation of two before 
the Raman process (cf. Fig. ||b). They are surrounded 
by a ring of atoms which originates from those site with 
an atomic occupation of one before the Raman transi- 
tion (see Fig. gc). Finally, Fig. gd shows the superfluid 
molecular density after ramping the optical lattice down. 

We now turn to the description of the coherent Ra- 
man transitions involved in creating molecules in step 
(iii) of our scheme. There are several constraints on the 
choice of detunings and laser intensities (see Fig. g) to 
maximize SI. First, the detuning A from an allowed ex- 
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FIG. 3: 2D lattice, a) Initial atomic SF density |$a| 2 . b) 
Number of atoms n a and c) number of molecules m immedi- 
ately after Raman conversion, d) Final SF molecular density 
$f,| 2 - The chemical potential is /i = 2. 5 Jo, the additional 
trap potential is given by ei = 2Jq(x\ + j/ 2 )/5a 2 , where Xi, %)i 
denote the coordinates of well i. The depth of the optical 
lattice is changed 10 times slower than in Fig. [j]. 



cited photoassociation resonance |^pr) (cf. Fig. ||) must 
be large compared to the natural linewidth jpr of the 
intermediate state to suppress spontaneous Raman scat- 
tering. Also, one should not detune half-way between 
two resonances since interference between these two in- 
termediate states will lead to a minimum in the effective 
two-photon Rabi frequency O. For appropriate choices of 
A and the intermediate state |\Ppr), the coupling to all 
intermediate states other than |\&pr) can be neglected. 
Second, as outlined above, we have to ensure that the 
process where a loi photon is absorbed before a u\ pho- 
ton (schematically shown in the inset of Fig. ^) has neg- 
ligible probability. This process causes trap excitations 
of single atoms and thus leads to heating. If both of 
these conditions are fulfilled the effective Rabi frequency 
SI on Raman resonance for the first Raman step is given 
by fi = $liQ2/2A, where SI1.2 are the Rabi frequencies 
for the first (second) step, and the effective spontaneous 
emission rate is 7 = 7prQi /4A 2 . 

For the case of 87 Rb, the following Raman pathway is 
viable for producing X l T^{v — 0, J = 0) molecules: 



X x E+(« tr 
X^+(v 



ap-0) 

= 120) - 



■* A 1 E u (v = 213) -» X r E+(v = 120) 
A 1 Yi u (v = 185) -> -X" x E+(w = 52) 
X^liv = 52) -► A 1 Y, u (v = 24) -» A" x S+(t; = 0). (3) 

A two-step Raman pathway also exists for producing 
a 3 E+(v = 0, J = 0) molecules. For the first step of 
the singlet pathway, the vibrational spacing near the 
A 1 E tl (u = 213) level is 110 GHz, lPR » 12 MHz, and the 
binding energy of the X l E+(v = 120) is 31.9 GHz. Given 
a trap frequency v of 1 MHz and intensities 1\ — lW/cm 2 
and I2 = 10~ 3 W/cm 2 for the first Raman step of Eq. (^), 
one obtains Sli = 0.71 MHz and Sl 2 = 3.7 MHz. For a 



red detuning of 200 linewidths we get Q = 1.1 kHz and 
and an effective spontaneous Raman scattering rate of 
7 = 5.5 Hz. For the process described in the preceding 
paragraph one finds that the undesirable loi — u)% process 
has a blue detuning of 29.5 GHz off the same intermedi- 
ate level and a Rabi frequency that is more than a factor 
of 100 below SI. Although the pathway for producing 
singlet molecules is partially optimized, no attempt to 
optimize the intensities and detunings has been made. 
However we have chosen values to show that the undesir- 
able u>2 — process can be sufficiently suppressed. The 
rate limiting step in the overall pathway is determined by 
the matrix element of the first Raman step. The matrix 
elements for all subsequent steps are at least 3-orders of 
magnitude larger than this one. 

For experimental purposes, molecules in higher lying 
vibrational levels of the ^ 1 E+ or a 3 E+ state should 
also allow the formation and detection of a molec- 
ular BEC. Vibrational relaxation caused by inelastic 
molecule-molecule collisions, that ultimately limits the 
lifetime of the resulting molecular condensate, will in gen- 
eral be strongly suppressed as the kinetic energy in the 
exit channels increases. This occurs as the vibrational 
spacing increases and should allow for observations of 
molecular BEC for moderately bound vibrational levels. 

In conclusion, generation of a MI phase of atoms al- 
lows efficient conversion of atoms to molecules, and to 
obtain a molecular condensate via "melting" in a quan- 
tum phase transition. This idea can be immediately gen- 
eralized to e.g. heteronuclear molecules, or one could use 
laser chemistry to build more complex composite objects 
(trimers etc.) and, possibly, corresponding condensates 
by quantum melting. The key to designing these pro- 
cesses is the fact that the MI phase provides us with a 
given small number of particles per lattice site (reaction 
partners) whose few body dynamics can be understood 
in all detail and controlled via laser interactions. 

Discussions with E. Ticsinga, P.S. Julienne and 
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in part by the Austrian Science Foundation, EU Net- 
works, and by the U.S. Office of Naval Research. 



[1] 
[2] 



[3] 



[4] 

[5] 
[6] 
[7] 
[8] 



J.R. Anglin and W. Ketterle, Nature, 416, 211 (2002). 
R. Wynar e t ai. Science 287, 1016 (2000); P. D. Drum- 



ond-mat/01 IO574 K - Goral, M. Gajda, and 



mond et al. 

K. Rzazewski, Phys. Rev. Lett. 86, 1397 (2001). 
M. Mackie, R. Kowalski, and J. Javanainen, 
Phys. Rev. Lett. 84, 3803 (2000); J. Javanainen 
and M. Mackie, Phys. Rev. Lett. 88, 90403 (2002). 
M. Holland, J. Park, and R. Walser, Phys. Rev. Lett. 86, 
1915 (2001). 

C. Gabbanini et al, Phys. Rev. Lett. 84, 2814 (2000), 

D. Jaksch, et ai, Phys. Rev. Lett. 81, 3108 (1998). 
M. Greiner, et al, Nature 415, 39 (2002). 

H.-J. Miesner et al., Science 279, 10051007 (1998). 



■5 



[9] For a review see: T. Calarco, et al. Fortschritte der 
Physik, 48, 945 (2000). 
[10] M.P.A. Fisher et al., Phys. Rev. B 40, 546 (1989); K. 
Sheshadri et al, Europhys. Lett. 22, 257 (1993); J.K. 
Freericks and H. Monien, Europhys. Lett. 26, 545 (1994); 



L. Amico and V. Penna, Phys. Rev. Lett. 80, 2189 (1998). 
[11] F. H. Mies, E. Tiesinga, and P. S. Julienne, Phys. Rev. 
A 61, 022721 (2000). 



